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Abstract. This paper presents a millimeter wavelength aperture-synthesis study of the spatial variations of the 
chemistry in the envelope around the deeply embedded low-mass protostar L483-mm on ~1000 AU (5") scales. 
Lines of 8 molecular species including CN, C ls O, CS, C 34 S, HCN, H 13 CN, HCO+ and N 2 H+ have been observed 
using the Owens Valley Radio Observatory Millimeter Array. Continuum emission at 2.7-3.4 millimeter is well-fit 
by an envelope model based on previously reported submillimeter continuum images down to the sensitivity of 
the interferometer without introducing a disk/compact source, in contrast to what is seen for other protostellar 
objects. A velocity gradient in dense material close to the central protostar is traced by HCN, CS and N2H + , 
and is perpendicular to the large-scale CO outflow, with a pattern consistent with rotation around a ~1 Mq 
central object. Velocity gradients in the propagation direction of the outflow suggest a clear interaction between 
the outflowing material and "quiescent" core. Significant differences are observed between the emission morpholo- 
gies of various molecular species. The C 18 interferometer observations are fit with a "drop" abundance profile 
where CO is frozen-out in a region of the envelope with temperatures lower than 40 K and densities higher 
than 1.5 xlO 5 cm -3 , which is also required to reproduce previously reported single-dish observations. The N2H + 
emission strongly resembles that of NH3 and is found to be absent toward the central continuum source. This 
is a direct consequence of the high CO abundances in the inner region as illustrated by a chemical model for 
the L483 envelope. The observed CN emission forms a spatial borderline between the outflowing and quiescent 
material probed by, respectively, HCO + and N2H + , and also shows intermediate velocities compared to these two 
species. A scenario is suggested in which CN is enhanced in the walls of an outflow cavity due to the impact of 
UV irradiation either from the central protostellar system or related to shocks caused by the outflow. 

Key words, individual objects: L483-mm, stars: formation, ISM: molecules, ISM: abundances, astrochemistry 



; 1. Introduction l2002t I2004J) these high-resolution observations make it 

• possible to resolve and address the spatial variation of 

The chemistry of star-forming regions shows a richness the cnemistry in thc pro tostellar envelope on -1000 AU 

and complexity reflecting large variations in the physical ^ scaleg _ 0nly a few protostars have previously been 

conditions found in these environments. For example, the studied in tMg chemical detail and the comb ination with 

thermal evolution of the pre- and protostellar cores results line rac i iat i ve transfer models and information about thc 

in evaporation and freeze-out of molecules illustrating the larger gcale structure from s i n gl e -dish observations makes 

important interplay between the solid-state and gas-phase a unique discussion about the chemical differentiation at 

chemistry. Recent single-dish surveys of low-mass prot 0| varying temperatures and densities possible. 
stars iUorgensen et all l2002t l2004d ISchoier et alJ l2002t 

iMaret et al l2004h have illustrated that the chemistry may Recent single-dish continuum studies of embedded pro- 
be severely affected by the thermal history of cores, with tostars have proved useful to establish the physical struc- 
radial variations of the temperature due to the heat- ture of the protostellar envelopes. A detailed description 
ing from the central protostar. This paper presents high- of the density and temperature variation is crucial for 
resolution millimeter wavelength aperture-synthesis obser- calculations of the molecular excitation and for the in- 
vations of a wide range of different molecules toward the terpretation of the chemical small-scale structure from 
embedded low-mass protostar L483-mm. In combination high resolution interferometer observations. Typically thc 
with a detailed radiative transfer model l|Je)rgensen et all envelope density and temperature structure can be con- 

strained from radiative transfer m odeling of subm illimc- 

Send offprint requests to: Jes K. j0rgensen ter continuum obs ervations fe.g.. Ijprgensen et "all 120021: 

Correspondence to: joergensen@strw.leidenuniv.nl ISm"rlev et alJ [20021 or from infrared extinction measure- 
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ments (e.g., lAlves et alJ l200lt IHarvev et al.l 12001^1 . The 

detailed profiles can subsequently be used as input for 
Monte Carlo radiative transfer modeling of single-dish 
line emission to constra i n the molecula r abu n dances 
Ce.g.. IBergin et all 120021: Ijorgensen eTaD l2002t l2004d: 
ISchoier et alJl2002|) . The models can also be used to pre- 
dict images of both continuum and line emission that can 
be directly compared to interferometer observations con- 
straining the envel ope structure down to scales of a few 
hund red AU (e.g., I.Iorgensen et alJ l2004at ISchoier et, alJ 
120041) . The caveat of the single-dish studies is their low 
spatial resolution of 10-15" or worse, which is why these 
observations predominantly are sensitive to material on 
larger scales (>2000 AU) unless high-excitation lines 
are observed. Recent interferometric studies of millime- 
ter continuum emission around low- mass protostars (e.g., 
Hogerheiide et alJll999t lLoonev et alJl200Ct IHarvev et all 



Table 1. Parameters for L483 from Ijgirgensen et alJ 
(12002ft . 



2OO3Uj0rgensen et alJl2004al ISchoier et alfl2004h have il 
lustrated the potential for probing the small-scale physical 
structure of the envelopes and constrain the presence of 
unresolved emission, possibly originating in circumstellar 
disks in these deeply embedded stages. 

Images of molecular line emission at similar resolu- 
tions can be used to discuss the detailed env e lope ch em- 
ical structure. For example. I.Torgensen et al.l 1 2004al) re- 
ported high-resolution observations of a range of molecular 
species toward the class protostar, NGC 1333-IRAS2A. 
The 3-6" (600-1200 AU) observations were interpreted 
in the context of envelope models constrained by single- 
dish continuum maps an d multi-transition line qbserva - 
tions from the survey by I.Torgensen et al.l ((2002L l2004ch . 
It was found that the single-dish envelope model could 
be successfully extrapolated to the smaller scales, lending 
further credibility to the approach and derived physical 
and chemical properties. 

L483-mm (IRAS 18148-0440; in the following just 
L483) is similar to NGC 1333-IRAS2A being a deeply 
embedded, low-mass protostar with a low bolometric tem- 
perature of w 50 K. In contrast to NGC 1333-IRAS2A, it 
shows a remarkable, asymmetric structure in the SCUBA 
maps. This asymmetry is likely to be the cause of the 
rather flat density distributions (n oc r~ p ) found in ra- 
diative transfer mode ling of the co ntinuum emiss i on by 
IJorgensen etaD (|2002ft (p = 0.9) and lShirlev et alJ l|2002j) 
(p = 1.2). The envelope parameters for L483 are given in 
Ta ble [J 

IPark et alJ l|2000h reported BIMA observations of 
HCO+ 1-0 and C 3 H 2 2 X2 - l i toward L483. They found 
two characteristic velocity gradients in the two species: an 
east-west velocity gradient in HCO + coincident with the 
larger scale CO outflow and a north-south velocity gradi- 
ent in C3H2 which they associated with global contraction 
of an envelope of a few thousand AU size. L483 has also 
been m apped in NH3 , most recently by Fuller & Wootter] 
(2000) who presented NH 3 maps from the VLA. They 
found that the overall asymmetry observed with SCUBA 
also shows up in the NH3 maps, but with NH3 lacking 
emission close to the protostar, which they ascribed to 



Distance, d 


200 pc 




9 Lq 


Thai 


50 K 


Envelope parameters: 


Inner radius (T = 250 K), B4 


9.9 AU 


Outer radius, -Riok 


l.OxlO 4 AU 


Density at 1000 AU, n(H 2 ) 


l.OxlO 6 cm -3 


Slope of density distribution, p 


0.9 


Mass, Miok 


4.4 Mq 



optical depth effects. They also noted a characteristic ve- 
locity pattern across the maps and suggested the presence 
of infall close to the central protostar. 

In addit i on to clear signatur e s in CO emission 
llFuller et all Il995l: lHatchell et ail llQQflt iTafalla et all 
120001) . an outflow driven by L483 is also seen through 
near-infrared emission, possibl y caused by scattered emis- 
sion in the outflow cavities teodapdll993 iFuller et alJ 
1199.1) . The absence of significant enhanceme nts of CH 3OH 
and S iO, as found in other outflows, led ITafalla et alJ 
(2000) to suggest that the L483 outflow is more evolved 
than those from other class objects. Since the driving 
source itself appears de e ply embedded ({Fuller et alJ ll995: 
IFuller fc Woottenll2000|) . ITafalla et all therefore suggested 
that L483 is in transition from the class to the class I 
stage. 

This paper presents high resolution 3 mm inter- 
ferometer observations of a wide range of molecular 
species toward L483. The main objective is to address 
whether the observed t rends in the single-dish survey by 
Ijprgensen et alJ l)2004cj) can be related to the variations in 
the chemistry in the L483 envelope. The paper is laid out 
as follows: Sect.|5|describes the details of the observations. 
Sect. 13 discusses the continuum maps of L483 and com- 
pares to the predictions from th e envelope mode l base d 
on the SCUBA observations from ljorgensen et al.l l)2002j) . 
Sect. 0] presents the line observations. Sect. brings to- 
gether the results for the line emission and discusses the 
implications for the chemistry in the environment of L483. 



2. Observations 

L483-mm (a 2 ooo = 18 h 17 m 29!8, c5 2 00 o = -04°39'38"3) 
was observed using the Owens Valley Radio Observatory 
(OVRO) Millimeter Array 1 from October 2001 to 
February 2003 in three settings including HCN, H 13 CN 
and HCO+ (3.4 mm), N 2 H+, CS and C 34 S (3.2 mm), 
and C ls O and CN (2.7 mm) as summarized in Table [21 
Each of the settings was observed in the C and E con- 
figurations providing projected baselines ranging from 3 
to 45 kA. Compared to the BIMA HCO + maps presented 



The Owens Valley Millimeter Array is operated by the 
California Institute of Technology under funding from the US 
National Science Foundation (grant no. AST-9981546). 
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bv IPark et al ] <HH|, the (u,v) coverage of the tracks in 
this paper includes longer baselines but not the shortest 
spacings. These observations are therefore less sensitive 
to extended emission but provide higher resolution. The 
correlator gives spectral resolutions of w 0.15-0.2 km s _1 
over 128 channels covering each line with RMS noise lev- 
els of 0.1 Jy beam -1 per channel. The complex gains were 
calibrated by regular observations of the nearby quasar, 
j 1743-038, approximately every 20-25 minutes. The band- 
pass was calibrated by observations of strong quasars, the 
fluxes using the same quasars as secondary calibrators and 
Uranus and Neptune as primary calibrators. Calibration 
of th e data was performe d using the MMA reduction pack- 
age l|Scoville et al.lll993[) . 

Table 2. Summary of the observed lines. 



Molecule 


Transition 


Frequency 


H ia CN 


l-0 a 


86.3402 


HON 


l-0 a 


88.6318 


HCO+ 


1-0 


89.1885 


N 2 H+ 


l-0 a 


93.1737 


C 34 S 


2-1 


96.4129 


CS 


2-1 


97.9810 


C 18 


1-0 


109.7822 


CN 


1-0° 


113.4910 



Notes: a Hyperfine splitting observed. 




60 40 20 -20 -40 -60 60 40 20 -20 -40 -60 
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Fig. 1. Continuum maps at A=3.4, 3.2 and 2.6 mm from 
the OVRO observations compared to the SCUBA 450 /mi 
map. The contours are in levels of 2 a as given in Table |3J 
The position of the radio source from lBeltran et aTl ([2001) 
is indicated by the black '+' in the upper left, A = 3.4 mm, 
map. 



3. Continuum emission 

The continuum maps clearly show the central protostar 
close to the phase center (Fig. . The emission appears 
resolved and roughly follows the morphology from the 
larger scale SCUBA maps. Results of Gaussian fits to the 
vi sibility data are given in Table |3J The envelope model 
of Ijgirgensen et all l|2002j) predicts total integrated fluxes 
ranging from 0.16 to 0.32 Jy at wavelengths from 3.4 to 
2.6 mm. Fig. [5] shows that the envelope model reproduces 
the observed fluxes on all the baselines, but the total flux 
recovered in the interferometer maps is only 5-10% of the 
total model flux. This is likely due to the lack of sensitivity 
in the interferometer observations to extended emission. 

From Fig. [5] it is also evident that very little emission is 
picked up on baselines longer than 20-25 kA (i.e., at scales 
less than 10"). Fig. [31 shows the fitted point source fluxes 
after subtracting the envelope model (F compac t, listed in 
Table • It can be seen that the emission is consistent 
with the zero-expectation level (the expected amplitude 
signal due to noise alone in the absence of source emis- 
sion) within la (3.2 and 3.4 mm observations) and 2a 
(2.6 mm) . This contrasts the si tuation for NGC 1333- 
IRAS2A llJoreensen et alJl2q04al) and IRAS 16293-2422 
and L1448-C l)Sch5ier*^t*alTl2004j) which show compact 
emission not explained by the envelope models. In fact the 
upper limit to the compact flux at 2.6 mm restricts the 
mass of a possible disk to < 0.04 M Q assuming a temper- 



ature of 30 K and optically thin emission. R adio measure- 
ment s of L483 at centimeter wavelengths ijBeltran et alJ 
l200l|) show a central source of 0.20 mJy (6 cm) and 
0.31 mJy (3.6 cm). These radio observations are compared 
to the millimeter data from this paper in the insert in 
Fig. |21 As can be seen, a spectral index consistent with 
optically thin emission cannot simultaneous explain the 
millimeter and centimeter observations; however, a flat, 
positive spectral index would be consistent with the emis- 
sion at centimeter wavelengths and the limits on the mil- 
limeter observations. This would be the case if the flux 
from centimeter to millimeter wav elengths is from therm al 
free-free emission as suggested bvlBeit ran et al.ll|200l|) . 

4. Line emission 

4.1. Morphology 

Maps of the total integrated emission of each of the ob- 
served lines are shown in Fig. 01 All species are detected to- 
ward the continuum position and most are stretched in the 
east-west direction of the elongated core from the SCUBA 
maps. 

The N2H + emission is seen to follow the contours 
of the dust emission from the SCUBA continuum maps 
closely, except toward the position of the central star 
where the emission appears suppressed. This gives N 2 H + 
emission a unique "peanut-shape" compared to the re- 



4 J.K. j0rgensen: Imaging chemical differentiation around the low-mass protostar L483-mm 

Table 3. Results from fitting Gaussians to the visibility data for the three continuum datasets. 



Wavelength [mm] 


3.4 


3.2 


2.6 


RMS [mJy beam" 1 ] 


0.5 


0.9 


0.8 


Beam size (HPBW) ["] 


7.1x6.5 


8.8x7.6 


6.2x5.2 


ftotai [mJy] 


13±2 


12±2 


11±2 


X-offset ["] 


-0.5±0.9 


0.1±0.8 


1.5±0.3 


Y-offset ["] 


-0.9±0.7 


-0.3±0.5 


-0.9±0.3 


^compact™ [ m Jy] 


1.0 ±0.5 (2.0) 


4.8 ± 1.0 (5.2) 


5.0 ±0.9 (4.5) 



"Fitted emission from a compact source after subtraction of the model prediction for the envelope emission. The numbers in 
parentheses indicate the zero-expectation level, i.e., the expected amplitude signal due to noise alone in the absence of any 
source. 
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Fig. 2. Plots of visibility amplitudes vs. projected baseline line length for the three continuum datasets at A=3.4, 3.2 
and 2.6 mm. The solid circles indicate the observations and the lines the predictions from the model for the L483 
envelope. 



maining species, with the strongest emission at two lobes, 
respectively east and southwest of the continuum peak 
position. C 18 in contrast is mainly seen very closely con- 
fined to the region around the central protostar. The emis- 
sion is resolved and interestingly shows an almost triangu- 
lar shape similar to that found in the SCUBA maps (e.g., 
lower right panel of Fig. ^| : the two shortest sides in the 
north-south direction on the western side of the core and 
in the east-west direction on the northern side of the core, 
and the longest side stretching from the northeast to the 
southwest. 

HCO + shows the most extended emission. It has a 
slightly curved structure around the central protostar, be- 
ing almost anti-correlated with the N2H+ emission. A sim- 
ilar "twist" can be seen in the single-dish C O maps of 
lHatchell et ail <|l999() and lTafalla et all (|2000|) and in gen- 
eral the HCO + appears closely correlated with the out- 
flow, as was also conclude d based on the BIMA HCO + 
maps bv IPark et ail l)2000i) . Fig. compares the 2MASS 
All-Sky Quicklook 2 K s image of L483 with the observed 

2 The Two Micron All Sky Survey (2MASS) is a joint project 
of the University of Massachusetts and the Infrared Processing 
and Analysis Center/California Institute of Technology, funded 
by the National Aeronautics and Space Administration and 
the National Science Foundation. The presented image was 
obtained from the NASA/IPAC Infrared Science Archive, 



HCO + emission. Near-perfect overlap is seen between the 
infrared and HCO + emission in the western outflow lobe. 
This in agreement with the suggestion by iFuller et all 
l|l995|) that the infrared emission is light from the cen- 
tral protostar scattered off the outflow cavity walls. The 
HCN and CS emission appear to be more centrally con- 
centrated around the central protostar but do show traces 
of the same elongated east-west feature. In addition a 
north-south feature perpendicular to the main elongation 
is seen in both species extending over « 20". Their weaker 
isotopic species, C 34 S and H 13 CN, only show unresolved 
emission around the central continuum peak. 

CN shows a very narrow feature elongated in the east- 
west direction. Toward both N 2 H + lobes the CN emission 
makes a twist as seen in the HCO + maps, but for CN the 
emission is much narrower. In fact CN is found to trace 
the boundary between the N2H + and HCO + emission (see 
Sect. 15.2} . The two CN hyperfme components otherwise 
show similar structures indicating quite homogeneous ex- 
citation conditions along the main elongation. 

As found for the continuum observations, the interfer- 
ometer picks up only a fraction of the line emission seen by 
the single-dish telescope. Fig. compares the single-dish 

which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology, under contract with the National 
Aeronautics and Space Administration. 
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40 20 -20 -40 40 20 -20 -40 40 20 -20 -40 
RA offset ["] RA offset ["] RA offset ["] 

Fig. 4. Overview of all observed lines at OVRO. The black line contours indicate the line emission integrated from 2.5 
to 8.0 km s _1 for CS, HCN, HCO + , and N2H + and from 4.0 to 6.5 km s _1 for the remaining species with levels given 
in steps of 2a. The grey line contours indicate the SCUBA 450 /im emission, the grey-scale map indicates the 3.4 mm 
continuum emission. For HCN, H 13 CN and N2H + the emission is integrated over the main hyperfine component. The 
CN*1 and CN*2 maps show the emission of the two CN 1-0 hyperfine components at 113.489 and 113.491 GHz, 
respectively. 



spectra from the Onsala 20 m telescope i|j0rgensen et alJ 
l2002ll2004rft with those from the interferometer datacubes 
convolved with Gaussians similar to the size of the single- 
dish beam. The interferometer spectra recover only a frac- 
tion of the emission close to the systemic velocity of the 
cloud (w 5.3 km s _1 ) whereas that in the wings agrees 
better, i.e., is less subject to resolving out. 

The large degree of resolving out, together with the 
asymmetry of the source, naturally complicates the in- 



terpretation of the line observations and, in particular, 
renders direct calculations of, e.g., column densities based 
on the absolute values of the interferometer observations 
impossible. The interferometer observations do, however, 
give the location and sizes of the brightest emission and 
one can thereby utilize comparisons between the maps of 
the different species to address the spatial variations of 
the chemistry. Also, predictions from the line radiative 
transfer model can be used for comparison to the inter- 
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10 15 



LSR 



Fig. 6. Comparison between the single-dish spectra (grey) and spectra extracted from the interferometer datacubes 
convolved with the single-dish beam (black). Each spectrum from the interferometer data has furthermore been 
multiplied by a factor 5 to make comparison between the lineshapes easier. As in Fig.0]the spectra marked CN*1 and 
CN*2 shows the emission of the two CN 1-0 hyperfine components at 113.489 and 113.491 GHz, respectively. 



ferometer observations, especially for species that show a 
relatively simple structure such as C 18 (see discussion in 
Sect. EH . 



4.2. Velocity field 

Velocity gradients are seen for a number of species to- 
ward L483 as illustrated in channel maps of the lines of 
HCO+, HCN, CS and N 2 H+ in Fig.0 The HCO+ maps 
clearly show the twist of the outflow toward the east for 
the red-shifted lobe, with distinct cores at offsets (15",- 
5") at 6.6 km s^ 1 and the edge of the channel maps at 
> 6.6 km s _1 . The first of these two clumps are also clearly 
seen in the CS channel maps at the same velocity. This 
core is merging with the central protostellar core at offsets 
(3", -3") seen toward all four species. In the maps of the 
blue-shifted emission the northern part of the central core 
is also seen at (0",5") for velocities around 4.2 km s _1 . 
The outflow emission is less obvious in the blue-shifted 
lobe. A last feature is a red-shifted core at (-15", 0") in 
HCN, HCO+ and CS at velocities of 6.3 km s" 1 . The two 
peaks of the N2H+ emission are seen to be unrelated to 



emission in the other maps and peak at respectively 5.1 
and 6.3 km s _1 . 

As discussed by I.Torgensen et al it may be 

somewhat problematic directly inferring velocity gradi- 
ents from interferometer maps with missing short-spacings 
since the interferometer is predominantly sensitive to ma- 
terial at the more extreme velocities. As can be seen from 
Fig. El the emission in the line wings is indeed least sub- 
ject to resolving out compared to the emission close to the 
systemic velocity, but in general the observed line profiles 
appear similar from both single-dish and interferometer 
observations indicating that they are still probing much 
the same material. 

For N2H + , the hyperfine splitting further compli- 
cates the interpretation of the velocity field. As was also 
found for NGC 1333-IRAS2A, however, each component 
is rather narrow, not affected by the outflow, and well- 
represented by single Gaussians. To derive the velocity 
field probed by N2H + , the emission from the 7 hyperfine 
components was therefore fit simultaneously to the spectra 
in each pixel. The lincwidths for the individual hyperfine 
components were taken to be constant and fit together 
with the systemic velocity field, the overall normalization 
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Fig. 7. Channel maps for the HCO + , HCN, CS and N2H" 1 " emission. Contours are given in steps of 3<t. 



of the entire hyperfme group and the relative intensities of 
the hyperfme components. The fitted velocity field is illus- 
trated in Fig. [S] An overall red-blue asymmetry is found 
between the eastern and western lobes. Toward the source 
position it is seen that the velocity gradient is directed 
toward the north-south as in the cases of HCN and CS. 
The ver y characteristic two lobes in N2H + are also seen in 
NH 3 bv lFuller fc Woottenl l|2000() . who furthermore found 
that the NH3 emission shows a two velocity component 
structure with the NE lobe red-shifted to 5.5-5.6 km s _1 



and the SW lobe blue-shifted to 5.3-5.4 km s" 1 . This 
is exactly what is found for our N2H + data, indicating 
that the N H3 and N2H" 1 " emis sion probe the same ma- 
terial. Both iFuller et all l|l995|) and iTafalla et alJ (|2000h 
proposed L483 as a core being broken up by the action 
of the outflow resulting, for example, in the bipolar near- 
infrared nebulosity. The velocity field seen in N2H" 1 " and 
NH3 and its similarities with the outflow velocity pattern 
would suggest that such an interaction in fact is ongoing 
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Fig. 3. Compact (unresolved) emission as a function of 
wavelength at 2.6, 3.2 and 3.4 mm. For each data point 
the expected flux in the absence of any point source is indi- 
cated by an open star. Fitted power-law distributions with 
spectral indexes of 2, 3 and 4 are shown with the dashed, 
dotted and solid lines, respectively. The insert show the 
millimeter fluxes co mpared to the centimeter data from 
iBeltran et alJ l|200lh . Here the dashed line indicates the 
flux from centimeter wavelength with a spectral index of 
0.8. The solid line indicates a spectral index of 3, corre- 
sponding to optically thin emission from dust with opaci- 
ties n v oc v$ where (3 = 1 at millimeter wavelengths. 
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Fig. 5. 2MASS near-infrared K s image of L483 (grey-scale 
and white contours) compared to the integrated HCO + 
1-0 emission (black contours; in steps of 3cr). Notice the 
near-perfect agreement between the infrared nebulosity 
and the outflow emission probed by HCO + . 
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Fig. 8. The velocity field from the N 2 H + maps as colored 
contours compared to the integrated emission (black line 
contours). The velocity has been derived by fitting the 7 
hyperfine components of the N2H 4 " 1-0 line over the entire 
map. 



although these species are not directly probing the out- 
flowing gas. 

iFuller k. Woottenl l)2000|) found the most blue-shifted 
emission toward the central star with velocities of 5.1- 
5.3 km s -1 , and suggested that the more blue-shifted 
emission toward the central star is a result of infall. In 
our N2H+ maps it is also seen that the emission has its 
most blue-shifted peak about 5" northwest of the pro- 
tostellar mm/cm source. However, an accompanying red- 
shifted component is present on the other side of the cen- 
tral source (i.e., to the southeast). This is further sup- 
ported by our HCN and CS maps which show a similar 
vel ocity pattern aro und the central protostar, also seen 
bv lPark et alJ l(2000h in C3H2 maps. A tempting sugges- 
tion for this velocity gradient is rotation around the cen- 
tral protostellar object and the propagation axis of the 
outflow. The CS and HCN lines have rather high criti- 
cal densities and appear to probe only the densest gas on 
small scales close to the central protostar. Likewise C3H2 
is typically only found in the dense envelope material not 
affected by the outflow (e.g. , iBachiller fc Perez Gutierrez! 
ll997tlTafalla fc Mverslll997l . 

Fig. El shows the position- velocity diagram for the CS 
emission. The coordinate system has been rotated to have 
the Y-axis in the direction of the north-south velocity gra- 
dient around the central protostar (i.e., 12° from north 
through east) and the X-axis in the direction of the out- 
flow perpendicular to this. The largest velocity gradient 
is seen in the north-south direction of the CO outflow. At 
velocities ~ 1 km s^ 1 from the systemic velocity the gra- 
dient is found to be consistent with linear expansion of 
0.27 km s _1 arcsec -1 or 2.6 xlO~ 4 km s _1 pc _1 at the 
distance (200 pc) of L483. 
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Fig. 9. Position-velocity diagram for the CS emission ro- 
tated 12°, so the Y-axis is in the direction of the velocity- 
gradient in the dense material close to the central proto- 
star (see also Fig. 1100. The dashed lines indicate Keplerian 
rotation around a 1.2 Mq central object in the upper panel 



and linear expansion of 0.27 km s 
panel. 



1 arcsec 1 in the lower 



A clear velocity gradient is also observed in the north- 
south direction: plotted on top of the position- velocity dia- 
gram is the predicted rotation curve for Keplerian rotation 
around a 1.2 M© central object. This velocity gradient 
fits the data at large velocities compared to the systemic 
velocity, but at velocities of 4.5-5 km s _1 the emission 
appears to stretch southwards in contrast to the predic- 
tion for Keplerian rotation. The cause of this is the offset 
in the north-south direction between the red and blue- 
shifted lobes of the outflow. Alternative explanations to 
the pure Keplerian rotation (e.g., a combination of rota- 
tion and infall such as sug gested for another embedded 
low-mass YSO, TMC1, bv lHogerheiideil200l|) cannot be 
ruled out based on the data presented in this paper, since 
the confusion with the outflow around L483 complicates 
the interpretation. 



5. Discussion 

5.1. Thermal structure, depletion of CO and resulting 
chemistry 

The most striking feature in the maps of the integrated 
emission in Fig. 0] is the characteristic "peanut-shaped" 
N 2 H + maps, with N2H + being absent from the gas-phase 
close to the central protostar. This likely reflects destruc- 
tion of N 2 H+ by reactions with CO, as su ggested to be 
the ca se in p rotostellar enviro n ments by , e .g..lBergin et al 
i|200ll) and Ijprgensen et all l)2004alch . Ij0rgensen et al 
( 2004cJ) found a clear anti-correlation between the abun- 
dances of CO and N 2 H + , which was suggested to be caused 
by CO depletion, and thereby less N 2 H + destruction, in 
the cold part of the envelope. It is therefore interesting to 
compare the N 2 H + and CO maps presented in this paper: 
as further illustrated in Fig.[IJ3l the N 2 H+ and C 18 peaks 
are clearly anti-correlated with CO being present towards 
the central protostar. This supports the claim that the 
C 18 maps peak toward the regions where CO has evap- 
orated from the grain mantles whereas N 2 H + peaks away 
from the central continuum position where CO is frozen- 
out. 

As discussed in Sect. 14.11 a significant fraction of the 
line emission is resolved out. Since the structure of the 
C 18 emission appears relatively simple, we can calcu- 
late the emission from the envelope ad opting the phys- 
ical st ructure and CO abundances from l.lorgensen et alJ 
(2002) as shown in Fig. ^2 This model is compared to 
"jump" models where the CO abundance increases at radii 
where the temperature is above a given threshold mimick- 
ing evaporation of CO. As can be seen from Fig. 1111 the 
constant abundance model provides a better fit to the in- 
terferometer data: the "jump" models either overproduce 
the observed emission at intermediate baselines (models 
with T cv = 20 K and 30 K) or do not produce a large in- 
crease at the sh ortest baselines (model with T f v = 40 K). 
As discussed in l.lorgensen et al ] ll2002t l2004r| . however, 
although the constant CO abundance and jump mod- 
els can successfully explain the intensities of the higher 
J = 2 — 1 and J = 3 — 2 CO lines, they generally under- 
produce the low J = 1 — emission. It is therefore not 
unexpected that such models have problems explaining 
the radial variation of the observed CO emission. 

As an alternative, I.Torgensen et alJ l(2004bO suggest 
a "drop" model where CO is frozen out only in regions 
where the temperature is lower than the evaporation tem- 
perature T cv , while at the same time the density is higher 
than a given density, ride, so that the depletion timescale 
is shorter than the age of the core. This is similar to 
the case seen for pre-stellar co res fc.g.. lCaselli et al .111 99<j 
iTafalla et al 1 l2002t iLee et al.l E)0^ where CO is frozen 
out towards the core center where the density high, and 
the corresponding timescale for depletion thereby shorter 
than the lifetime of the core. Such "drop" models can ex- 
plain both high and low J CO lines, and the difference in 
abundances between more and less massive envelopes seen 



10 



J.K. j0rgensen: Imaging chemical differentiation around the low-mass protostar L483-mm 




40 30 -20 

RA offset ["] 




40 20 -20 -40 

x-offset ["] 



Fig. 10. Comparison between N2H 4 " vs. C 18 emission. 
Upper panel: The integrated N2H" 1 " and C 18 emission 
(black and white contours) plotted over SCUBA 450 /im 
image (grey-scale). The flattened direction of the core (i.e., 
12° from the east-west axis) has been indicated by the 
dashed line. The coordinate system with the rotation ve- 
locity gradient (Fig. 0) has been indicated by the dotted 
lines. Lower panel: the N2H" 1 " and C 18 emission in a strip 
in the flattened direction of the core (dashed line in upper 
panel). Note the clear anti-correlation between the peaks 
in C 18 and N 2 H+. 

in ljdrgensen et aTl l)2002j) . and they work very well in de- 
scribing the radial v ariations of molecules such as H2CO 
l|Schoier et al.ll2004h . 

Models for the single-dish CO emission were calculated 
for L483 and compared to the line intensities reported in 
Ijorgensen et all ll2002l) . A drop model with n^c = 1.5 x 
10 5 cm -3 and T ev = 40 K, an undcpleted CO abundance 
of 2.7xl0 -4 and depletion in the drop region by a factor 
50 was found to provide good fits to the observed single- 
dish emission as indicated by the reduced % 2 for the fit 




( u 2 + v Z)l/2 [ kX ] 



Fig. 11. Comparisons between the observed visibilities 
(solid circles) and model predictions (grey lines) for the 
C 18 emission averaged over the observed line (4.0 to 
6.5 km s _1 ) plotted vs. projected baseline length. Model 
1-3 are jump models with T cv of 20, 30 and 40 K, respec- 
tively. Model 4 is a drop model with T cv = 40 K and 
ride = 1-5 x 10 5 cm -3 . Model 5 is a constant abundance 
model. Note that only model 4 i s consistent with all sing le- 
dish line intensities presented in ljoreensen et aD l|2002h as 
shown in Tabled 

as given in Table 0| This model has been compared to 
the interferometer data in Fig. ^2 an d is found to give a 
good fit to the observed C 18 1-0 emission at practically 
all baselines. Freeze-out up to a temperature of 40 K is 
needed in order not to produce too much emission at the 
intermediate baselines and the inclusion of an outer low 
density region, where the CO abundance increases, nicely 
reproduces the observed trend at the shortest baselines. 

5.1.1. Chemical implications 

The CO abundance structure has important chemical im- 
plications. As described above, distinct zones with and 
without depletion of CO also account for the character- 
istic shape of the N?H + emission. In their NH3 maps 
iFuller & Woottenl ( 2000) found a very similar structure 
with a "valley" devoid of emission in the north-south di- 
rection through the position of the protostar but otherwise 
following the larg e scale dust continu um emission from 
the SCUBA maps. ICaselli et aP l(2002i) likewise suggested 
a close correlation between emission seen in N2H 4 " and 
NH3 based on a large single-dish survey of N2H + emis- 
sion toward dense cores. The precursor for both N2H + 
and NH3 is N2 which due to its low binding energy does 
not deplete unless at very high ages ari d densities (e.g. 
iBergin fc Langerll997tlBergin et al.l2002|) . NH3 maintains 
a balance with the overall chemical network of N-bearing 
species, including N2H + , predominantly through reactions 
with NH^, and at low (depleted) CO abundances these 
species establish a closed network. At "standard" CO 
abundances, reactions between C + and NH3, however, 
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Model 






Xrcd,s 


Xrcd.i 


1: 


Jump T*ev 


= 20 K 


221 


478 


2: 


Jump T ev 


= 30 K 


18.9 


8.0 


3: 


Jump T cv 


= 40 K 


10.8 


4.7 


4: 


Drop T cv 


= 40 K; ride = 1.5xl0 5 cm" 3 


1.5 


1.3 


5: 


Constant 




8.0 


3.1 



a Reduced \ 2 for fits to single-dish line intensities. ^Reduced x 2 for fits to interferometer data as shown in Fig. 1111 



drive the nitrogen into H2NC + and from there to CN, 
HNC and HCN, whereas the N2H + is destroyed through 
reactions with CO. 

These trends are illustrated with a chemical model for 
the L483 envelope using the chemical net work and ap- 
proach described in lDotv et alJ I2O O2I 12004 |) and ado pting 
the physical structure from Ijgirgensen et alJ l)2002j) with 
CO depleted in the zone constrained by the single-dish and 
interferometer observations above. As shown in Fig. 1121 
this gives exactly an increase in N2H + and NH3 over the 
region where the CO is depleted. Note that although the 
CO abundance goes up again in the outermost part of 
the envelope and the NH3 and N2H+ abundances con- 
sequently drop, this occurs at distances of « 8000 AU 
corresponding to 40" from the compact source in Fig. 
which is close to the boundary of the N2H + emission. 
This is, however, probably not the main reason for the 
outer edge to the N2H + emission: it is more likely caused 
by the critical density of the N2H + 1-0 transition which 
is w 10 5 cm -3 , a density which is also reached at this 
distance from the central source. The N2H + emission is 
therefore not sensitive to the outer region with low den- 
sity material where CO is undepleted. 

Of course a detailed radiative transfer model for the 
N2H + and NH3 emission could address some of these is- 
sues in more detail, such as it was done for the CO emis- 
sion. The problem is, however, that the spherical model 
clearly breaks down for the scales probed by the N2H + 
emission. In fact the asymmetry may well be the reason 
that N 2 H + is seen in so clearly distinct zones away from 
the central protostar in the L483 case, since N 2 H + be- 
comes marginally optically thick toward the source center 
in an entirely spherical core. Therefore N2H + will not be 
sensitive to the innermost region where the CO comes off 
the dust grains in the symmetric case. 

5.2. UV irradiation of outflow cavity walls 

HCO + is found to be most prominent towards the out- 
flow as is typically seen at large and small scales (e.g., 
iHoeerheiide et all Il998|) . whereas the N 2 H+ emission 
probes the cold part of the quiescent cloud. As illustrated 
in Fig. If 31 the CN emission appears in a boundary region 
between HCO + and N2H+, i.e., between the outflowing 
and quiescent cloud material. A possible explanation for 
this "borderline property" of the CN emission could be 
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Fig. 12. Molecular abundances as function of radius 
for L483, adopting the chemic al model appropr iate for 
IRAS 16293-2422 at 10 4 years l|Dotv et alJl2004|) . Upper 
panel: the density and temperature profile for L483 shown 
with dashed and solid lines, respectively. Lower panel: the 
abundances of CO, NH 3 , N 2 H+ and HCN from the model 
calculations. The CO abundance (solid line) structure is 
assumed to follow a "drop" abundance structure with an 
evaporation temperature, T cv , of 40 K and a depletion 
density, n^c, of 1.5xl0 5 cm -3 as indicated by the two ver- 
tical arrows. The N2H+ (dashed line) and NH3 (dotted 
line) abundances are found to increase by one to two or- 
ders of magnitude over the region where CO is depleted, 
whereas HCN (dashed-dotted line) is seen to follow the 
CO abundance, dropping by a factor 5-10 over the same 
region. 

that CN is enhanced in the walls of an outflow cavity 
pr obed by the HCO + emiss ion. 

JOrgcnscn et al found a clear trend between 

the CN abundances and CS/SO abundance ratio and sug- 
gested that this might reflect the strength of ultraviolet 
radiation in the protostellar envelope. L 483 is the source 
in the sample of lJorgensen et a l. (2004c) with the highest 
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Fig. 13. Comparison between HCO+, N 2 H+ and CN 
emission. Upper panel: SCUBA 450/im image with HCO + 
emission and N 2 H + emission as white and black solid 
line contours, respectively, and the vertical centroid of the 
CN emission overplotted (black/white dashed line). Lower 
panel: emission of the three species in the vertical direction 
averaged over the eastern outflow lobe (offsets between 10 
and 25" in the coordinate system aligned with the flat- 
tened direction of the core (dashed lines in Fig. 1171)1 . 



ratio of CN abundances in the inner region relative to the 
outer region of the envelope (i.e., [CN] 32 /[CN]io, where 
[CN]32 and [CN]io are the CN abundances inferred from 
the high excitation 3-2 lines at 340.24 GHz and low ex- 
citation 1-0 lines at 113.49 GHz, respectively). For most 
sources in the sample this ratio is less than 1 indicating 
that the CN abundances are enhanced predominantly in 
the low density material in the outermost region of the en- 
velope. For L483, the ratio is 2.5 indicating a higher CN 
abundance in the denser material. If the CN abundances 
are probing the strength of the UV field, this argues in 



favor of an internal source of the UV radiation in the case 
of L483. 

Such correlations between the CN abundances and UV 
field have been suggested to be prese nt in gaseous disks 



aro und more evolved class I /II objects fchitrev et al. 
pi l200li iBergin et all l2003t Ivan Zadelhoff et alJ 
Thi et al.l 120041) which are known to be strong emitters 



of UV radiation resulting from the ongoing accretion. In 
the deeply embedded stages, however, the UV radiation 
cannot escape far out into the envelope due to the large 
extinctions. This will prevent direct detection of any UV 
flux from the central star-disk system, but this emission 
may still be dominant at small scales close to the cen- 
tral protostar. It is therefore of high interest to search for 
specific chemical probes that may be used to address the 
impact of UV in the inner envelope and thereby poten- 
tially also the ongoing accretion. 

The UV radiation field may furthermo re be of large 
impor tance in cavities cleared by outflows. ISnaans et alJ 
(|l995|h for example, suggested that the presence of high 
J — 6 — 5 CO lines toward embedded protostars is re- 
lated to heating of material in the outflow cavity walls 
by the UV radiation escaping from the central proto- 
star. The_jie^£jJifrared nebulosity toward L483 (Fig. [3] 
and iFuller et alJ l|l995|n suggests that such a cavity is 
present. Alternatively the enhanced UV radiation could 
be a result of the shocks associated with the L48 3 outflow 
which are directly prob e d by the H2 emission l|HodaniJ 
ll994HBuckle et al.lll994 - lMolinari et alJ l|200lh suggested 
that recombination in the post-shocked gas could illumi- 
nate outflow cavities behind Herbig-Haro objects produc- 
ing a PDR resulting: in diffuse C + emission seen by ISO 
l|Molinari et alJl200lt iMolinari & Noriega-Crelrooll2002h . 

To get a handle on this possibility a simple estimate 
of the extent of the CN emitting region can be made: 
the CN emission is unresolved in its transverse direction, 
meaning that the extent of the emissio n must be < 5". 
In the model for the L483 envelope from |j0rgensen et alJ 
l|2002() the density at 20" from the protostar is n(H 2 ) = 
3 x 10 5 cm~ 3 . Introducing the extent of the CN emission 
as I < 5", an upper limit to the column density of the CN 
emitting region can be estimated: 



V H , < 4 x 10 21 cm" 2 



n(H 2 ) 



3 x 10 5 c 



(1) 



or a visual extinction, Ay < 4. Models o f pho- 
ton dominated reg i ons (e.g., Ijansen et al.1 Il995i 
ISternberg k. Dalgarnol Il995|) indicate that this is the 
extinction range over which CN would be enhanced pre- 
dominantly due to photodissociation of HCN. Alternative 
explanations such as excitation effects appear less likely, 
since CN is the only molecule showing this peculiar 
morphology. 

If CN in fact is tracing the interface region between 
the outflow and the quiescent core probed by HCO + and 
N 2 H+, respectively, this should also be reflected in the 
observed velocity field of the three species. As discussed 
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Fig. 14. Velocity field in HCO+, N 2 H+ and CN emission: 
average centroid velocity of the three species in strips per- 
pendicular to the outflow axis across the eastern part of 
the core. 

above the velocity field probed by N2H + shows a red- 
blue spatial asymmetry aligned with that seen in HCO + , 
suggesting that the "quiescent" material probed by N2H+ 
is being accelerated by the outflow. Fig. El compares the 
velocity field in the eastern lobe of the three species along 
the flattened direction of the core (dashed line in Fig. I10|) . 
The HCO + does indeed show the largest velocity gradient 
across the region with less steep increases in CN and N 2 H + 
toward the lower densities in the core. Again CN is seen 
to be intermediate between the outflowing and quiescent 
material in support of the suggested scenario. 

6. Conclusions 

We have presented an analysis of the physical and chem- 
ical properties of the class protostar L483-mm on 500- 
1000 AU scales. This paper demonstrates how the ID en- 
velope model derived from single-dish continuum and line 
data can be applied to the interpretation of the millimeter 
interferometer observations. The main conclusions are as 
follows: 

1. The continuum e mission at 3 mm is well- fit by the en- 
velope model from lj0rgensen et alJ (|2002l) . In contrast 
to the case for other recently studied class objects 
(e.g., NGC 1333-IRAS2A, L1 4 48-C and IRA S 1629 3- 
2422; I.Tctrgensen etaD j2004a|) : ISchoier et all l|2004l) L 
an additional point source (e.g., a disk) is not needed 
to fit the observed continuum emission. Assuming an 
optically thin 30 K disk, this puts an upper limit to 
the disk mass of 0.04 M Q . 

2. The C ls O emission is foun d to be well-described in 
a "drop abundance" model l|j0rgensen et al.ll2004blcfl 
where the abundance is high in the outermost regions 
(densities lower than 1.5xl0 5 cm~ 3 ) and innermost 
(temperatures higher than 40 K) parts of the envelope 
but frozen out in between. 



3. N2H+ and C 18 are found to be clearly anti- 
correlated, with C 18 centered on the central contin- 
uum source whereas N2H + has two distinct peaks away 
from the central star. This is interpreted as the com- 
bined effects of CO freeze-out at low temperatures and 
destruction of N2H + by reactions with CO. The N2H + 
emission resembles that o f NH3 , previously reported by 
iFuller fc Woottenl l|2000|) . which can be understood if 
NH3 is destroyed through reactions with C + when the 
CO abundance is high, and as clearly illustrated by a 
chemical model for the L483 envelope. 

4. HCN and CS emission probe the dense material close 
to the central protostar. A velocity gradient perpen- 
dicular to the outflow propagation directio n, also rec- 
ognize d in previous C3H2 observations bv IPark et alJ 
( 200c| , is interpreted as rotation around 1 Mq central 
protostar. 

5. CN is found to trace a boundary between the quies- 
cent material probed by N2H + and the larger scale 
outflow seen in HCO + . A possible explanation is that 
CN probes material in the outflow cavity walls (seen 
as an infrared nebula in 2MASS K s images) where 
its abundance is enhanced as a result of UV irradi- 
ation. This is further supported by the observed veloc- 
ity field: the gradient introduced by the outflow is seen 
in both HCO+, CN and N 2 H+, with HCO+ showing 
the largest velocities relative to the systemic velocity 
followed by CN (in the cavity walls) and N2H+ (in the 
core material). This also suggests that a clear inter- 
action between the outflowing and quiescent material 
around L483 is taking place. Possibly the outflow is in 
the process of dispersing the protostellar core, causing 
its characteristic asymmetric shape. 

This paper illustrates the potential of high-resolution 
millimeter interferometer observations for addressing the 
spatial variations in the chemistry around protostellar ob- 
jects. In particular, future high resolution observations 
from the SMA, CARMA and eventually ALMA will con- 
firm or reject the drop abundance model through imaging 
of the radial structure of the emission from high excita- 
tion lines of, e.g., C 18 0, probing more uniquely the higher 
temperatures and densities in the envelope. Also high res- 
olution observations of high excitation lines of CN will 
make it possible to further address the importance of UV 
radiation in the inner envelope close to the central pro- 
tostar and confirm the relation between CN abundances 
and strength of the UV field. In this context direct imag- 
ing of lines of, e.g., atomic carbon in the high frequency 
windows will also be important to understand the specific 
chemical effects. Such observations will become feasible 
with the SMA and ALMA. An additional test would be 
deep mid-infrared Spitzer observations: for example, deep 
imaging of the thermal emission in the outflow cavities and 
searches for established probes of the UV field such as the 
presence and properties of specific emission lines and PAH 
features, provide complementary information which could 
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be used to establish a more detailed 2D model of complex 
protostellar sources such as L483. 
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